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Abstract 

We give a microscopic explanation of both Debye and non-Debye thermalization 
processes that have been recently reported by Gall and Kutner [1]. Due to reduc- 
tion of the problem to first passage phenomena we argue that relaxation functions 
f(t) introduced by the authors directly correspond to survival probabilities S(t) of 
particles in the considered systems. We show that in the case of broken ergodicity 
(i.e. in the case of mirror collisions) the survival probability decays as a power law 
S{t) = r/t. 
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1 Introduction: motivation and aims 



A few decades ago, it was widely accepted that the relaxation of a typical 
macroscopic observable A towards its equilibrium value is described by expo- 
nential function 



fit) = M£ = e-^ 



where t represents time and AA(t) = A{t) — A(oo). On that score, probably 
the best-known example is the Newton's law of cooling where the relaxation 
function f(t) refers to thermalization process (i.e. A = T). Recently, however, 
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it was discovered that a number of decay phenomena observed in nature obeys 
slower non-exponential relaxation (c.f [2,3,4,5,6]) 




(2) 



It was also noticed that non-Debye relaxation is usually observed in systems 
which violate the ergodicity condition. In such systems, the lack of ergodic- 
ity results from long-range interactions, microscopic memory effects or from 
(mult i) fractal structure of phase space. It was argued that such systems are 
well described by the so-called nonextensive statistics introduced by Tsallis 
[7,8]. The conjectured relaxation function for such systems [9,10] is given by 
the g-exponential decay 1 



which is equivalent to the formula (1) for q ~ 1, whereas for q > 1 it coincides 
with (2). 

In the context of the ongoing discussion on possible relations between noner- 
godicity and nonextensivity, the phenomenon of non-Debye relaxation that has 
been recently reported by Gall and Kutner [1] seems to be particularly inter- 
esting (see also [13]). The authors have numerically studied a simple molecular 
model as a basis of irreversible heat transfer trough a diathermic partition. 
The partition has separated two parts of box containing ideal point parti- 
cles (i.e. ideal gases) that have communicated only through this partition (see 
Fig. la). The energy transfer between the left and right-hand side gas samples 
has consisted in equipartition of kinetic energy of all outgoing particles collid- 
ing with the partition during a given time period. The authors have analysed 
and compared two essentially different cases of the system's dynamics: 

i. the first case, where the border walls of the box and the diathermic parti- 
tion have randomized the direction of the motion of rebounding particles, 
and 

ii. the case, where mirror collisions of particles with the border walls and 
the partition have been considered. 

They have found that although the mechanism of heat transfer has been anal- 
ogous in both cases the long-time behaviour of both thermalization processes 



1 As a matter of fact, the only satisfactory proof of the property (3) exists for such 
systems in which nonextensivity arises from intrinsic fluctuations of some parame- 
ters describing the system's dynamics [11,12]. 




(3) 
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has been completely different. In the first case (i.) ordinary Debye relaxation 
of the system towards its equilibrium state has been observed 



AT(t) ~ e 



-t/r 



(4) 



where AT(t) = Ti(t) — T 2 (t) is the temperature difference between both gas 
samples, while in the second case (ii.) the power-law decay has been noticed 



In order to describe the phenomenon of the non-Debye relaxation Gall and 
Kutner [1] have derived an extended version of the thermodynamic Fourier- 
Onsager theory [14,15] where heat conductivity was assumed to be time- 
dependent quantity. The authors have argued that from the microscopic point 
of view the non-Debye relaxation results from the fact that the gas particles 
always move along fixed orbits in the case (ii.). They have also argued that 
this regular motion may be considered as nonergodic, violating the molecular 
chaos hypothesis (Boltzmann, 1872). 

In this paper we propose a more rigorous microscopic explanation of both 
Debye and non-Debye thermalization processes reported by Gall and Kutner. 



2 Microscopic model for non-Debye heat transfer 

In the paper [1], the authors have analysed two-dimensional systems consisting 
of two gas samples of comparable size (see Fig. la). Here, due to analytical 
simplicity we assume that one gas sample is significantly larger and denser than 
the second one i.e. the larger sample may be referred to as a heat reservoir 
with constant temperature = const (see Fig. lb). We also assume that 
the smaller sample is confined in a square box of linear size / and its initial 
temperature equals T . It is natural to expect that thanks to the existence of 
the diathermic partition the temperature of the smaller sample will tend to 
the reservoir temperature T(t) — > in the course of time. 

At the moment, before we analytically justify the relaxation functions (4) and 
(5) let us recall crucial assumptions of the numerical experiment performed 
by Gall and Kutner [1]. First, the authors have defined the temperature of 
the given gas sample T(t) as proportional to the average kinetic energy of all 
particles in the sample 



AT(t) 



T 



(5) 



t 



(6) 
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Fig. 1. (a) Original experimental system [1]. Two gas samples exchanging heat 
through a diathermic partition, (b) Reduced model system. Gas sample in thermal 
contact with a huge heat reservoir at constant temperature (detailed description 
is given in the text.) 

Second, they have assumed a monoenergetic energy distribution function P(e) = 
5(e — kTi) as the initial condition for each gas sample i — 1, 2 (see Fig. la). In 
fact, since the applied thermalization mechanism evens out kinetic energies of 
all particles colliding with the diathermic partition at a given time not only 
initial but also final (i.e. equilibrium) energy distributions are monoenergetic. 

Now, having in mind the above assumptions, one can simply conclude that 
during the heat transfer occurring in the system presented at Fig. lb the 
particles of the smaller and thinner gas sample get the final energy immediately 
after the first collision with the diathermic partition 2 i.e. 6q — > e^, where 
Eq = kT and = kToo (6). It is possible due to the existence of the huge 
and dense heat reservoir which causes that the number of particles with energy 
colliding with the diathermic partition at a given time is overwhelmingly 
larger than the number of particles with the energy e which at the same time 
collide with the partition from the other site (see the inset in Fig. lb). 

The above considerations allow us to write the temperature difference between 
both gas samples in the following way 

AT(t) = T{t) -T 00 = ^P(T - T^), (7) 



2 Gall and Kutner have proved that the systems presented at Fig. la possess a 
similar feature for asymptotic times (cf. Eq. (45) in [1]). 
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Fig. 2. Ideal point particles in the box with mirror border walls (detailed description 
is given in the text). 

where N(t) is the number of particles of the smaller sample which have not 
hit the diathermic partition by time t. Now, one can see that the relaxation 
function f(t) of the considered systems is equivalent to the survival probability 
S(t) = N(t)/N 

fit) = T ® ~ T °° = S(t). (8) 

J-O ~ -too 



The last formula makes us possible to reduce the phenomena of Debye and 
non-Debye relaxations to the first passage processes [16,17,18]. In this sense, 
the case (i.) of rough border walls directly corresponds to the problem of 
diffusing particles in a finite domain with an absorbing boundary. The survival 
probability S(t) typically decays exponentially with time for such systems [16]. 
That is the reason way the thermalization process characterizing the case (i.) 
is equivalent to Debye relaxation (see Eq. (4)). In the next paragraph we show 
that the case (ii.), where mirror bouncing walls and absorbing diathermic 
partition are taken into account, is indeed characterized by the power law 
decay of the survival probability S(t) (see Eq. (5)). 

In order to achieve the claimed scale-free decay of the survival probability S(t) 
one has to find the so-called first passage probability F(t) (i.e. the probability 
that a particle of the considered gas sample hits the diathermic partition for 
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the first time at the time t) 

t 

S(t) = 1 - J F(t')dt'. 



(9) 



At the moment, let us remind that before the collision with the diathermic 
partition each particle has the same velocity vo (i.e. mvl/2 = Sq = kT ), thus 
the distribution F(t) can be simply calculated from the particle path length 
distribution F(r) i.e. 



F(t) = F{r) 



dr 
~dt 



F(v t)v , 



(10) 



where r = v^t. Now, due to the symmetry of the considered problem (that is 
due to the equivalence of paths and — > A — > b — > 

c — >■ d, see Fig. 2a) one can deduce the following relation 

F(r) = l -P{a) 



da 
dr 



(11) 



where 0<a,/5<7r/2 and respectively sin a = x/r whereas sin (3 = (21 — x)/r 
(see Fig. 2b). The last formula expresses the fact that particles can either move 
to the left (i.e. towards the diathermic partition) or to the right (i.e. towards 
mirror reflection of the partition). Next, assuming uniform initial conditions 
P(a) = P(f3) = 2/ it in the long time limit (i.e. for small angles when sin a ~ a 
and sin (5 ~ 0) one obtains 

Fit) = — r 2 . (12) 

TTVq 



Finally, using the relation (9) one gets the desired power-law decay of the 
survival probability which justifies the non-Debye thermalization process (5) 

r 2/ 
S(t) = -, where r = . (13) 



We have numerically verified the last relation for a few different values of 
both the initial velocity Vq and the box size /. In all the considered cases we 
have obtained very good agreement of recorded survival probabilities with the 
formula (13) (see Fig. 3). 
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Fig. 3. Survival probability S(t) against time t in systems presented at Fig. 2. 
Points correspond to results of numerical simulations whereas solid lines represent 
theoretical prediction of the formula (13). 



3 Summary and concluding remarks 



In this paper we have given a microscopic explanation of Debye and non-Debye 
thermalization processes that have been recently reported by Gall and Kutner 
[1]. The authors have studied a simple molecular mechanism of heat transfer 
between two comparable gas samples. Owing to analytical simplicity we have 
reduced the problem to one gas sample being in thermal contact with the 
huge and dense heat reservoir at constant temperature. For the case we have 
shown that the thermalization mechanism described by Gall and Kutner can 
be reduced to first passage phenomena. Taking advantage of the idea we have 
found an analytical justification for both exponential (4) and non-exponential 
(5) relaxation functions observed in [1]. 
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